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Abstract 

We consider the process of photon decay in quantum electrodynamics with a CPT- 
violating Chern-Simons-Iike term added to the action. For a simplified model with 
only the quadratic Maxwell and Chern-Simons-like terms and the quartic Euler- 
Heisenberg term, we obtain a nonvanishing probability for the decay of a particular 
photon state into three others. 
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1 Introduction 



The propagation of light in Maxwell theory with an extra CPT-violating 
Lorentz-noninvariant Chern-Simons-like term in the action has been studied 
both classically [1] and quantum mechanically [2]. The Maxwell term and 
the Chern-Simons-like term are quadratic in the photon field. These terms 
combined lead to birefringence, even in empty space. Additional (Lorentz- 
invariant) higher-order photonic terms in the action could perhaps produce 
other effects such as the decay of photons. (See also Ref. [3] for a discussion 
of the decay of a massive Dirac fermion in theories with spontaneous Lorentz 
and CPT violation.) 
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The purpose of this paper is to study the process of photon decay in detail, 
expanding on the brief remark in Section 6 of Ref. [2] . In order to concentrate 
on the essentials, we keep the theory as simple as possible. In fact, the the- 
ory we start from is just quantum electrodynamics — the theory of photons, 
electrons and positrons (see, e.g., Refs. [4,5]). It is, of course, known that the 
quantum effects of the electron-positron field lead to a quartic coupling of 
the photons (giving, e.g., light-by- light scattering), which is described by the 
quartic part of the so-called Euler-Heisenberg Lagrangian [6] . 

The photon model considered in this article has only the Maxwell and quar- 
tic Euler-Heisenberg terms, together with a hypothetical spacelike Chern- 
Simons-like term which breaks Lorentz and CPT invariance. For the moment, 
the precise origin of this Chern-Simons-like term can be left open, but at least 
one possible mechanism has been identified [7,8] (see Ref. [9] for a review). 
A more extensive discussion of the possible origin and consequences of the 
Chern-Simons-like term can be found in Ref. [2] , which also contains further 
references. 

The outline of our paper is as follows. The model is presented in Section 
2, together with some basic facts on the polarization states of the "photons." 
The Lorentz noninvariance of the model allows for photon decay and the rel- 
evant kinematics is discussed briefly in Section 3. (Some technical details are 
relegated to Appendix A). The matrix element for a particular decay chan- 
nel is then calculated in Section 4 and the corresponding partial decay width 
for a photon at rest is given in Section 5. (Numerical results for the phase 
space integral are presented in Appendix B). For this particular case, there 
is, in principle, one other decay channel available, but it does not contribute 
as shown in Section 6 with details relegated to Appendix C. The total de- 
cay rate of the particular photon state at rest is discussed in Section 7. The 
corresponding mean life time is expected to be large but finite. 

2 Model 

The Lagrangian density of the model considered in this paper is 



where the free part Cq consists of the usual Maxwell term [4,5] and an addi- 
tional Chern-Simons-like term [1,2], 



Our conventions are {g^v) — diag(l, —1, —1, —1) and eoi23 = 1, together with 

n = c = i. 

The Chern-Simons-like term in Eq. (2.2) has a mass parameter m. (Exper- 
imentally, there are tight constraints [1,10] on this mass, m ^ 10~^^ eV, as 



(2.1) 




(2.2) 
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will be discussed further in Section 7.) The Chern-Simons-like term contains, 
in addition, a purely spatial "four- vector" 



(77'^) = (0,77^77^77^) , 



(2.3) 



in terms of a "three- vector" fj of unit length, I?)! ^ = 1. The parameters r/^ are 
fixed once and for all, hence the quotation marks. In fact, the condition of hav- 
ing fixed parameters (coupling constants) r]^ in the Lagrangian Cq effectively 
selects a class of preferred inertial frames; cf. Refs. [1,2,3]. 

In the following, we will also write the Maxwell field strength F^j, = S^^i/ 
— d^A^ in terms of the electric and magnetic fields, E^. = F^q and B^. = 
ekimFim/'^ with ei23 = 1. If the gauge needs to be fixed, we will use the radiation 
(Coulomb) gauge [4,5,11], 



d-A^O. 



(2.4) 



Further details on the photon propagation and the microcausality of the free 
theory (2.2)-(2.4) can be found in Ref. [2]. 

The interaction term Ci in Eq. (2.1) is taken to be the quartic Euler- 
Heisenberg Lagrangian [6] 



2a^ 
45 M4 

45 M4 



^\ 2 



{\E\^-\B\^y + 7 (E-B) 



(2.5) 



where M is a mass-scale and a dimensionless coupling constant. For quan- 
tum electrodynamics with electron mass M ^ 511 keV and fine-structure con- 
stant a ~ 1/137, precisely this term appears in the one- loop effective gauge 
field action (the corresponding Feynman diagrams have a single electron loop 
with four external photon lines). The quartic Euler-Heisenberg term in the 
effective action is relevant for photon energies much less than M, which will 
be the case for the process discussed in the present paper (with photon ener- 
gies less or equal to m). For further details on this effective action, see, e.g., 
Refs. [11,12,13] and references therein. 

The free Lagrangian Co in the radiation gauge ^ives rise to two modes of 
plane- wave solutions with polarization vectors e±{k) and dispersion relations 



u!±{k) — (^k^ + - ± - m V w? -\- Ak"^ cos^ 9 ^ 



1/2 



(2.6) 



with k =\k\ and kcosO = k ■ fj. For generic wave vectors /c, the unit vectors 



^1 = {fj — cos9 k)/ sm6 , ^2 = {k x rj) / sin 9 , k = k/\k\, 



(2.7) 



form an orthonormal tripod with positive orientation, ^ ^2 = k. The polar- 
ization vectors of the two modes are then given by 
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e±{k) 



2 uj± cos 9 ii ^ i{T ± m) ^2 



with 



r = r{k) = Vm2 + 4A;2 cos2 9 . 



(2i 



(2.9) 



As usual, the free quantized photon field can be defined in terms of creation 
and annihilation operators, 



1 



+ 



a+(k) e+{k) exp (^—icu+t + ik ■ + H.c. 



a^{k) e-(k) exp (^—iu-t + ik ■ +}i.c. ^ . (2. 



10) 



The creation and annihilation operators have standard commutation relations, 

[as{k),at{l)]^ [al{k),al{l)]^0, [a,{k),al{l)] ^ S,,t S\k - I) , (2.11) 

for s,t G {+,—}. The Fock space of photon states is then readily con- 
structed [11,12,13], starting from the vacuum state |0) with the property that 
as{k) |0) = for any (fc, s). 

The free electric and magnetic field operators can be expressed in the same 
way as Eq. (2.10), replacing the polarization vectors e± there by 



f±(k) = iuj± e± + m (^e± ■ {k x fj)^ k — 



2r(r±m) 



X 



(2 c;| cos 9 ii T iLJ± (r ± m) 6 T m sin ^ (r ± m) k) (2.12) 



for E and by 

b±{k) = ik X e± — 



ik 



V2r(r ± 



m) 



(2uj±cos9 i2±i{r±m)ii) (2.13) 



for B. 

Before continuing with the quantized theory, it may be useful to consider the 
polarizations that result from Eqs. (2.12) and (2.13) for the classical electric 
and magnetic fields of a plane wave propagating in vacuo (see also Ref. [9]). 

For 9 = (i.e., a wave vector k in the preferred direction 17), the + and — 
modes are circularly polarized. But the vectors ^1 and ^2 from Eq. (2.7) are not 
well-defined at ^ = 0. Instead, they may be defined by a limiting procedure, 
where the limit k± I is performed for a wave vector k = k±ej_ + k\\fj . Here, 
e± is an arbitrary unit vector perpendicular to fj. The resulting vectors are 
^1 = — e_L and ^2 = e± x r). Consequently, the + mode is right-handed {R, 
negative hehcity) and the — mode is left-handed (L, positive helicity). 
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For < ^ < /t/2, both modes are elliptically polarized and the hehcities 
remain the same {R/L for +/— )• As long as kcos6 ^ m, the modes are 
essentially circularly polarized. In the opposite limit k cos 9 ^ m, both modes 
approach linear polarizations. 

For 9 — 7r/2, the two modes behave somewhat differently. The — mode 
becomes exactly linearly polarized, with /_ ~ iuj-^i ± k. But the electric field 
for the + mode has also a longitudinal component, /+ ~ "^^+^2 + n^-k. This 
longitudinal component is, however, relatively unimportant for frequencies 
UJ+ S> m. 

For 7r/2 < 9 < both modes are again elliptically polarized, but the 
helicities have changed {L/R for +/— ). For 9 = tt, finally, both modes are 
circularly polarized, with helicities opposite to the = case. 

We now return to the quantized theory. The energy-momentum "tensor" of 
the free theory is given by 

1 TT) 

^ _ gi^-pp'^Fp, - F^^F-'-g,, - - rft^'^'^^ApF,, . (2.14) 

This tensor is conserved, d^Q^" = 0, as may be checked with the help of the 
equations of motion, 

m 

d^F^'' ^ - e-^P^-rj^F^r . (2.15) 

The energy-momentum tensor (2.14) is, however, not symmetric, as long as m 
is nonzero. 

For the case of a purely spatial Chern-Simons "four- vector" rj^, the corre- 
sponding field energy and momentum are 

S = Jd^xe^^'^ld^x^^lEl^+lBl^) , (2.16) 



V' = J d^xQ^' = J d^x (e X b)' + ^fi' (b ■ a) 



(2.17) 



These operators in terms of the quantized free photon field must, of course, 
be normal-ordered [11,12,13]. One finds after some algebra 

:£: ^ j d^k [oj^{k) a\{k)a^{k) + uj^{k) a}_{k)a^{k)\ , (2.18) 
■.V: ^ j d^k [k a\{k)a+{k) + k a}_{k)a-{k)\ . (2.19) 

The energy eigenvalues, in particular, are nonnegative for the functions uj± as 
given by Eq. (2.6). 

Note, finally, that the Lagrangian (2.2), with fixed parameters = (0,?}), 
is translation invariant but not rotation invariant. This implies that energy- 
momentum is conserved but not angular momentum. 
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3 Kinematics of photon decay 



In this section, we study the kinematics for the decay of one initial photon 
with momentum q and energy u}±{q) into three final photons with energies 
ijj±{ki), i — 1,2, 3, and total momentum A;i+A;2+A;3 = q. The theory considered 
has been given in Section 2. In a first reading, it is possible to skip ahead to 
the last paragraph of this section, which summarizes the results. 

For the usual Lorentz-invariant case (m = 0), the total energy of the final 
state is minimal when all final momenta fcj are parallel to the initial momentum 
q. The reason is that adding perpendicular momenta, ki ki + A/cj with 
Aki-q — and J2i Aki — 0, increases the final energy, since ijj{ki-\-/S.ki) > u{ki). 

The same does not hold in our case (m ^ 0). For infinitesimal Ak with 
Ak ■ k — 0, we find 

Lu±{k + Ak) ^ uj±{k) ± (k ■ f/) (Afc • r)) , (3.1) 

which can be larger or smaller than uj±{k). This may then result in the lowering 
of the total energy of a given 3-photon final state, which can be shown as 
follows. 

Assume that we start from a configuration where the final momenta A'j, 

i = 1,2, 3, are parallel to the initial momentum q, which has q ■ fj ^ and 
q X i] ^ 0. Take the final particles 1 and 2 to have the same polarization (+ 
or — ) and leave the polarization of the final particle 3 unspecified (±). Then 

— * — * 

add infinitesimal perpendicular momenta Aki to the final momenta ki and 
choose Aki = —Ak2 and AA;3 = 0, so that overall momentum conservation 
is maintained. According to Eq. (3.1), one of the first two final energies is 
increased and the other decreased (the third final energy is, of course, un- 
changed). Moreover, the total final energy for the ( ±) case is decreased if 

the energy corresponding to the smaller momentum is decreased (for example, 

— * — * — * — * 

ii \ki\ < I/C2I then |Aa;_(/ci)| > |Aa;_(A;2)| ). The same result is obtained for 
the (+ + ±) case, as long as \kj ■ fj\ > m/\/2 for j = 1, 2. 

Still, it is worthwhile to study the case of parallel (collinear) final momenta 
in some detail. The frequencies u± may then be treated as functions of scalar 
variables, u± = ui±{ki) for ki = \ki\. Assuming both ki and k2 to be nonzero, 
one finds for the — state that 



{ki + ki) > UJ- {ki) + (/ca) , (3.2) 

— * — * 

except for the special case of /ci • 77 = /c2 • ''7 = 0, for which Eq. (3.2) becomes 

an equality. As a consequence, the decay > is allowed kinematically. 

This also implies that the decay H ^ is allowed, because a;+(g) > 0J-{q). 

For the + state, on the other hand, the following inequality holds: 

a;+(A;i + A;2) < a;+(/ci) + a;+(/c2) , (3.3) 



6 



again assuming nonzero ki and k2- In this case, no immediate conclusions can 
be drawn from the inequahty, because of relation (3.1). The inequalities (3.2) 
and (3.3) are proven in Appendix A. 

Another kinematically allowed decay is H > -\ . This can be shown by 

establishing the following inequality for appropriate parallel momenta: 

uj+{ki + k2 + fcs) > ^^+(^1) + (jJ-{k2) + uJ-ik^) , (3.4) 

with the notation ki = \ki\. Unlike the inequalities (3.2) and (3.3), inequality 
(3.4) does not hold for arbitrary values of the k^. Instead, some restrictions 
must be imposed on the momenta. 

Concretely, assume that /c2 -C m, A;2 -C ki and k^ <^ m, /ca -C /ci. Expanding 
the left-hand side of Eq. (3.4) up to first order in k2 + ks and the right-hand 
side up to first order in k2 and k^, we find 

uj^^{ki) ki + = ) {k2 + ks) > sm 9 {k2 + k^) . 

\ \l rn? + Akl cos^ 

Defining x = m/ki, this inequality can be written as 

2a;cos^6' ^/ x"^ x /—^ : 

1 + = > sin 6' IH h - Vx^ + 4 cos^ 9 

Vx^ + 4cos^e \ 2 2 J 

Both sides of the last inequality are manifestly positive. Squaring both sides 
and re- arranging them somewhat, we arrive at 

n o ^ 8x^ cos^ 9 x^ — 4 cos^ 9 

-x^ + 2 cos^ 9 + — — > X 



x^ + 4:Cos'^9 V^TIco^' 

This final inequality certainly holds if the left-hand side is larger than zero 
and the right-hand side smaller than zero. Writing 6 = 4 cos^ 9, the condition 
for the left-hand side gives 



x'^ <b{b- l)/4 + (6/4) Vb^ -2b + 9 
and the condition for the right-hand side 
x''<b. 

Both conditions can be fulfilled simultaneously, as long as cos^ 9^0. We 
might, for example, restrict the momenta, so that 6 > 2, which corresponds 
to cos^^ > 1/2. Then both conditions are equivalent and give x < \/2, which 
corresponds to ki > m/ y/2. [For the special case of cos^ ^ = 0, the quantity uj+ 
(c<j_) is effectively the energy of a massive (massless) particle, which implies 
that the inequality (3.4) does not hold for small enough values of ki/m.] 

The above discussion shows that the energy inequality (3.4) holds for certain 
restricted values of the final coUinear momenta. By continuity, this energy 
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inequality holds also if the final momenta deviate slightly from coUinearity. 
Therefore, our restrictions are restrictions to certain regions in phase space 
and not restrictions to submanifolds of measure zero. As a consequence, the 

decay H > -\ is allowed kinematically. 

Altogether, we have shown that the three decay channels — — > , 

H > and H — > -\ are allowed kinematically. For the five other decay 

channels, we have not been able to find allowed regions in phase space, either 
analytically or numerically. These allowed regions in phase space perhaps do 
not exist, but this remains to be proven. 



4 Matrix element H > 



We now calculate the matrix element for the decay of a + polarization 
mode into three — modes. To lowest order, there arc two contributions to this 
matrix element from the quartic Euler-Heisenberg Lagrangian. The first term 
in Eq. (2.5) gives the following contribution: 

{0\a+{q) : J d^x [\E{x)\^ - \B{x)\^y : al{ki)al{k2)aUh) \0) = 
5[uj+{q) - uj^{ki) - uj^ih) - uj^{ks) ] S^{q- ki - k2 - h) 



-1/2 



uj+{q) uj^{ki) uj^{k2) uj^{k3 

Uh) ■ Mh) - Uki) ■ b-{k2)\ [Mh) ■ niq) - b-{h) ■ bXiq) 

+ [Uh) ■ f4h) -b-{ki) -b-ih)] [f-{k2)-fl{q)-b.{k2)-bl{q) 

+ • f-{h) - Uh) ■ Uh)] [f-{ki) ■ mo) - Uh) ■ Ki^ 



(4.1) 



The second term in Eq. (2.5) gives, with different combinatorics, 

{0\a+{q) : J d'^x 7 {e{x) ■ B{x)f : al{ki)al{k2)al{h) |0) = 
1 



27r^ 



5[u+{q) - u^{ki) - u^{k2) - ou-iks) ] 5 (g - ki - k2 - ks) 



uj+{q) uj_{ki) uJ-{k2) uJ-ik^) 



-1/2 



X 7/4 



Uh) ■ L(^2) +/-(^2) • Hki)] [fl{q) ■ L(4) + /-(4) • fc;(g1 

+ [Uh)-Hh)+f-{h)-b-{ki)\ [fm-b-{k2)+f-{k2)-KiQ) 
+ [Uh) ■ b4h) + Uh) ■ b4k2)] [fm ■ b-iki) + Uki) ■ 



(4.2) 
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Unfortunately, the expressions involving the polarization vectors in Eqs. 
(4.1)-(4.2) arc rather complicated and one would like to simplify them by 
using some type of small m expansion. But such a procedure appears to be 
quite difficult and, in this paper, we keep the full expressions. 

The probabihty for decay of a single + state into three — states can be 
calculated by integrating the square of the amplitude over the final momenta. 
The result will, however, be a function of |g| and q ■ fj, for given m and M. 
For this reason, we turn to a simpler problem in the next section, namely the 
decay of a + state at rest. 



5 Pctrtial decay width for H — ^ 

The decay probabihty for + — > will be evaluated for the case of 

vanishing initial momentum, q — 0. The energy a;_|_(^ for g = equals m and 
the initial "photon" behaves like a massive particle at rest. The partial decay 
width is then obtained by squaring the matrix element of the previous section 
and integrating over the final momenta ki, for i = 1 ... 3. A combinatorial 
factor 1/3! = 1/6 must be inserted because of the three identical particles in 
the final state. 

Purely on dimensional grounds, we can write the partial decay width to 
order ct^ as 



-(4) 



512 TT^ 



K 



'2^ 
45 



m 

M8 ' 



(5.1) 



with a single number k > to be determined. 

The relevant phase space integral is defined as follows: 



I = j d^ki d^k2d^k3 6^{ki + k2 + k3)S[uj^{ki)+uj_{k2)+(jJ-{k3)-m] 



X 



a;+(0) uj-{ki) U-{k2) u;-{kz) g{ki, k2, k^) , 



(5.2) 



where the nonnegative function g depends on the electric and magnetic polar- 
ization vectors. Specifically, the factor g is given by the absolute value squared 
of the terms in large brackets in Eqs. (4.1) and (4.2) for ^=0, 



g{ki,k2,k3) 



/_(A;i) • f4k2) - 6-(^i) • ^'-(^2) /-(^3) • /;(0) + •••) + 4 



X 



/_(fci) • b.{k2) + f-ih) ■ b.{ki) /;(0) • b.iks) + 



(5.3) 



The expressions for /^(O) and &+(0) will be given in the next section. Note 
that only the mass-scale m appears in the integral (5.2), so that / oc m^. 
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A numerical calculation shows that the integrand of Eq. (5.2) is strictly 
positive over a finite region of phase space (see Appendix B). For the numerical 
constant k in the partial decay width (5.1), this implies 




which is the main result of the present paper. 



6 Pcirtial decay width for H — > -\ 

The decay of an initial + polarization state is in general more complicated 
than the decay of an initial — state. But for initial momentum q = a. ma- 
jor simplification occurs: the second decay channel + — > H no longer 

contributes, as will be demonstrated in this section. 

The matrix element for H > H decay follows from Eqs. (4.1) and (4.2) 

by replacing all — labels that occur in conjunction with ki by + labels. The 
energy delta function, for example, becomes 

5[u;+{q) - u;+{ki) - a;_(^2) - uj-ih) ] ■ (6.1) 

The energy conservation condition for the decay of a + photon at rest 
{q — 0) now takes the form 

(u+ih) - m) + cu-ih) + cu-ih) = . (6.2) 

— * 

This relation holds only if all three final momenta ki are zero, because each of 
the three terms in the above equation is positive semi-definite and zero only for 

— * 

vanishing momentum. The conditions ki — 0, for i — 1,2, 3, are restrictions to 
a lower-dimensional submanifold in phase space and have, therefore, measure 
zero. The resulting phase space integrals will lead to a zero contribution to the 
decay width, unless certain singularities of the phase- space integrand appear in 
the limit A;j — > 0. In the following, we will show that such infrared singularities 
are absent. 

The phase-space integral in spherical polar coordinates ki,9i,ipi, for i — 
1,2,3, takes the form 

/Q — * — * — * — * — * — * 
S {ki + k2 + ks) 5[u+{ki) +uj^{k2) + u^iks) - m] h{ki,9i,ipi) 

A;| /c| sin 62 sin 9^ d^ki dk2 dk^ d92 d9z d(p2 d(p3 
cu+(0) uj+{ki) uj_{k2) uj-{k-i) 

where the nonnegative function h{ki,9i,ipi) depends on the electric and mag- 
netic polarization vectors f±{ki) and b±{ki). Using the bounds uj^{k) > m and 
u}-{k) > ksm9 in the denominator and integrating over the momentum ki, 
we find for the phase-space integral 
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/"—*—* — * — * 

< 7+ — < / 5[<jj+{k2 + fca) + ^^-{k2) + uj-{k3) - m] h{kj, 9j, (pj) 

k2 k-s dk2 dk-i 662 dO^ d(p2 dcp^ 

where h now depends on the polarization vectors f±{kj) and b±{kj), for j ~ 
2, 3, and the vectors f±{—k2 — k^) and b±{—k2 — k^). 

The integral in Eq. (6.4) has two potential sources of infrared singularities, 
the first of which is the factor h{kj, 9j, (pj). But it is not difficult to see that 
h{kj,9j,(fj) is nonsingular in the limit kj ^ 0. As h depends only on the 
electric and magnetic polarization vectors, it suffices to demonstrate that these 
vectors are nonsingular in the infrared. Writing k = k±e± + k\\ f) with e± an 
arbitrary unit vector orthogonal to r), it can indeed be shown that 



and 



lini f+(k) — im(e± — ie± x 77) (6.5) 



lim f4k) = lim b+{k) = lim b^k) = . (6.6) 

fc— >0 fc— >0 fc— >0 



The second potential source of infrared singularities in the integral (6.4) 
is the energy delta function itself, if singularities of the type / dk6{k'^) are 
produced. The proof that this does not happen is somewhat involved and is 
relegated to Appendix C. 

Altogether, we find that the H > -\ channel does not contribute to + 

decay at vanishing initial momentum, ^= 0, so that 

ri'l+__ = o. (6.7) 

For g = 0, no other channels contribute, besides the two channels already 
considered. The reason is energy conservation (cf. Section 2): the energy of an 
initial + state at rest is m, whereas the energies of the final states + H — and 
+ + + are at least 2 m and 3 m, respectively. 



7 Total decay rate of a + photon at rest 

In this paper, we have studied the decay of photons in a relatively simple 
model, for which the Lagrangian (2.1) contains only the usual Maxwell and 
Euler-Heisenberg terms, together with a CPT-violating Chern-Simons-like 
term. The two mass parameters of the model are m, which enters the Chern- 
Simons-like term (2.2) linearly, and M, whose fourth power enters the Euler- 
Heisenberg term (2.5) inversely. The Euler-Heisenberg term has, in addition, 
an overall coupling constant a^. 

The photon of this model has two polarization states, labeled ±, with en- 
ergies and momenta given by Eqs. (2.6), (2.18) and (2.19). The total decay 
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rate to order «^ of a + state at rest follows from the sum of the partial decay 
widths (5.1) and (6.7), 



r^'^-^—K (^] — (71) 

with a constant k > 0. A numerical calculation shows that k is nonvanishing 
and, most likely, not vastly different from 1. [See Eq. (5.4) and Appendix B.] 
As explained in Section 2, the model considered is part of the effective gauge 
field action of quantum electrodynamics with an additional Chern-Simons-like 
term. An order of magnitude estimate for the photon mean lifetime can then 
be obtained by writing 



where for a and M the values of the fine-structure constant and the electron 
mass have been inserted and for m the experimental upper limit from polar- 
ization measurements on distant radio galaxies (see Refs. [1,10] and references 
therein). For m ~ 10~^^ eV as might be expected from the CPT anomaly [7], 
the photon lifetime would be larger by a factor of 10^^. 

The very large photon lifetime as indicated by Eq. (7.2) is perhaps not 
of direct relevance, at least for the current epoch in the history of the uni- 
verse. There remain, however, fundamental questions about the formulation of 
this particular CPT- violating extension of quantum electrodynamics if certain 
photon states are no longer absolutely stable (cf. Ref. [14]). 
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A Energy inequalities 

In this appendix, we prove the energy inequalities (3.2) and (3.3) for parallel 
momenta ki \ \ [The uj_ case has the additional conditions kj ■ fj ^ 0, for 
j — 1,2.] It is, in fact, not difficult to verify that these inequalities hold for 

— * 

sufficiently small kj = \kj\ <^ m. For the general case, we give a proof by 
contradiction. 

Assume that the inequalities do not hold for all values of ki and k2- Then 
there must exist values for ki and k2 with 

uj{ki + k2)^uj{ki)+cu{k2). (A.l) 
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[These special values of ki and /c2 could, of course, be different for the + and 
the — case. But the proof is analogous in both cases and we treat both cases 
at once, writing uj for either a;_|_ or a;_.] The equality (A.l) can also be written 
as 

k k 
Vphih + k2) = - — ^ Vpi,{ki) + - — ^ Vph{k2) , (A.2) 

fcl + K2 Kl + K2 

where fph(A;) = uj{k)/k is the absolute value of the phase velocity. 

Equation (A.2), now, holds only if one of the following three conditions is 
met: 

(i) Vpk{ki) = Vpi,{k2) = Vpk{ki + k2) , 

(ii) Vp^{ki) > Vpi,{ki + k2) and Vpi,{k2) < Vpi,{ki + k2) , 

(iii) Vp^{ki) < Vph{ki + /cg) and Vpi,{k2) > Vph{ki + fcs) • (A. 3) 

Each of these conditions implies that Vp\i{k) has an extremum somewhere in 
the momentum interval [min(A;i, k2), ki + k2]. But this conclusion contradicts 
the simple observation that, for both a;+ and a;_, the derivative of the phase 
velocity dvp^/dk is nonzero for all positive values of k. [The observation re- 

— * 

quires for the U- case the additional condition k ■ fj ^ 0.] This then implies 
that the assumption above Eq. (A.l) must be incorrect and that the inequal- 
ities (3.2) and (3.3) hold for all positive values of ki and /i;2, as stated in the 
main text. 

B Numerical result for / 

In this appendix, we report on a numerical calculation of the integral I , 

as defined by Eq. (5.2). A dimensionless quantity / is obtained by setting 

/ = / /m^. We can be relatively brief in describing our results since the 

calculation of a decay rate is well-known (see, e.g.. Sections 3.6 and 3.7 of 
Ref. [15]). We proceed in four steps. 

First, the integral over k^ is performed and the remaining six integration 
variables are taken to be the following dimensionless spherical coordinates: 

— * — * 

ki = \k-i\ , k2 = \k2\ , 9-1, 02, (p± = ip2±<Pl: (B.l) 

with polar angles 6'i, 6*2 G [0, tt] defined with respect to an axis in the preferred 
direction fj and azimuthal angles (fii,(p2 £ [0,27r] in the plane orthogonal to 
this axis. The angles </?+ and </?_ can be taken to run over [0, Air] and [— tt, -I-tt], 
respectively. 

Second, the energy delta function effectively sets ip- to a fixed value, 

(/?_ = (1 — 2n) arccosx, (B.2) 

for a known function % = x(A;i, 6*1, A;2, 6*2) and integer n = or 1. Also, there 
are the constraints that both \x\ and a;_(/ci,6'i) + a;_(/c2, 6*2) must be less 
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than 1, which can be implemented by introducing appropriate step functions 
G(l — Ixl) and 0[ 1 — uj-{ki, 9i) — cj_(/i;2, 62) ] into the integrand. 
Third, the resulting integral has the following structure: 

/= r dkidk2 r dOidOi t"" V X| ^ . tit, (B.3) 

Jq Jo 2 I ¥'-=(l-2n) circcosx;fe3=-fci-fe2 ' ^ 

with the integrand 

T=[A;i/u;_(A;i,^i)] [ A;2/a;_(A;2, ^2) ] e[ 1 - u;_(A;i, ^1) - u;_(A;2, ^2) ] 

X e (1 - I x(A:i, ^1,^:2,^2) I) \ sin ip_\-^ g{h,k2,h), (B.4) 

where and g are dimensionless functions (i.e., the expressions of the main 
text with m = 1). The integral (B.3) is complicated, but its integrand is still 
nonnegative; cf. the definition (5.3). 

Fourth, a numerical calculation with Mathematica [16] shows the inte- 
grand of Eq. (B.3) to be independent of and n. This effectively reduces the 
integral to a four-dimensional one and a numerical estimate gives 

I = I /m^ 0.2. (B.5) 

From this estimate, one obtains the result (5.4) quoted in the main text. 

The value (B.5) is to be considered preliminary. More work is needed to 
obtain an accurate result, both analytically (e.g., to make the independence 
of the azimuthal coordinate manifest) and numerically (e.g., to sample 
phase space efficiently) . 



C Analytic result for 7+ 

In this appendix, we demonstrate the vanishing of the integral /_|_ as 

defined by Eq. (6.3). We start with two preliminary steps. First, we introduce 
the following representation for the energy delta function: 

S(uj) = lim — ^ exp (-cu'^/a^) , a > . (C.l) 

Second, we replace the nonnegative function h{kj,6j,(pj) in the integral of 
Eq. (6.4) by the following bound: 

h{kj,e„ipj)< X: H{l,n)klk^, (C.2) 

l,n=0 

where H{l,n) are appropriate nonnegative numbers. This bound may be un- 
derstood from the observation that the factor h{kj,9j,ipj) is the absolute 
square of a sum of terms of the type 
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— * — * — * — * 

[vi{k2) ■ V2{k3)] [v3{-k2 - h) ■ M^)] > 



(C.3) 



where each vector Vi, for i = 1, ... ,4, is either an electric polarization vector 

— * — * 

f± or a magnetic polarization vector b±, as defined by Eqs. (2.12) and (2.13). 
Using the energy bounds u!+{kj) < kj+m and u!-{kj) < kj, each vector Vi{kj) 
can be shown to obey the following inequality: 



Vi{kj 



< 



vo,i{kj) + kjVi,i{kj) 



(C.4) 



where t?o,i and Vi^i are uniformly bounded vectors (i.e., they never exceed a 
certain length). In the product (C.3), then, the momenta k2 and ks show up at 
most quadratically, so that no powers higher than four occiir in the absolute 
square of a sum of such terms. This explains the bound (C.2). 
With these two steps, the bound (6.4) becomes 

poo rn r27r ^ 

</+__< lim — / dk2dk3 / (^^3 / d(p2d(p3 Y H(l,n) 
a-*o ay/Trm^ Jo Jo Jo 

X /C3+-^ exp^- a;+(A;2 + 4) - + ^-{k2) + ui-iks) ^ /a^^ . 

(C.5) 

Now change the momentum variables kj to kj/ a and use the estimates a;+ > m 
and [u}-{k2) +a;_(A;3)]^ > ou^{k2) +ou^{k3). This gives 



< /+__ < lim Y 



4 ^3+«+n 



dk2 dk3 dQ2 dB3 d^p2 d^p3 H(l, n) 4"^^ k^+^ 



X 



ew{-kl - ^ + ^ ^l + Acos^e2k^2/ml 



-^1 - ^ + ^ + 4 cos2 03 kl/ml ) , 
with rua = m/a. Since cos^Oj < 1, we finally arrive at 



(C.6) 



m 



l+l+n 



:,n=o 



3+/+n 



/ d02 d93 I d(p2 d(p3 Ii In , 
Jo Jo 

(C.7) 



in terms of the dimensionless integrals 



: J^^ dk k'+^ exp -HA + H^ ^1 + Ak^/ml ) . 



(C.8) 



The momentum integrals for / = 0, ...,4, can be evaluated analyti- 
cally with the help of Mathematica [16], but we are only interested in their 
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asymptotic behavior for rua — > oo. Making the change of variables y — k'^/ma 
and Taylor expanding the square root in Eq. (C.8), we obtain 

7^~Qm^'/^ for m„ ^ oo , (C.9) 

with positive coefficients q. 

Inserting the asymptotic results for Ii and /„ into the bound (C.7), we get 
(A is a positive constant) 

4 

< /+__ < A lim V H{1, n) m'+'+" q c„ m^'-'/'-"/' = , (C.IO) 

771 a, — ^OO 

l,n=0 

SO that 

(C.ll) 

which implies Eq. (6.7) in the main text. 
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